AFRL-APOSR-VA-IR- 2015-0232 


Geometric  Staictuie-Preseiving  Disc  letization  Sc  hemes  for  Nonlinear  Elasticity 


Atash  Yavari 

GEORGIA lECH  RESEARCH  CORPORAHON 


Oa/ 13/ 2015 
Final  Repoit 


DISIRIBUnON  A:  Distribution  approved  for  public  lelease. 


Air  Force  Research  Laboiatory 
AF  Office  Of  Scientific  Research  (APOSI^/  RIA 
Ariington,  Virginia  22203 
Air  Force  Materiei  Command 


FORM  SF  298 


Page  1  of 2 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing 
data  sources,  gathering  and  maintaining  the  data  needed,  and  compieting  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or 
any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Executive  Services,  Directorate  (0704-0188). 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information 
if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  fDD-MM-yyyy;  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  Toj 

14-08-2015  Final  Performance  01-06-2012  to  31-05-2015 

4.  TITLE  AND  SUBTITLE 

Geometric  Structure-Preserving  Discretization  Schemes  for  Nonlinear  Elasticity 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

FA9550- 12- 1-0290 

5c.  PROGRAM  ELEMENT  NUMBER 

61102F 

6.  AUTHOR(S) 

Arash  Yavari 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S]  AND  ADDRESS(ES) 

GEORGIA  TECH  RESEARCH  CORPORATION 

505  10TH  ST  NW 

ATLANTA,  GA  30318-5775  US 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES] 

AF  Office  of  Scientific  Research 

875  N.  Randolph  St.  Room  31 12 

Arlington,  VA  22203 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

AFRL/AFOSR  RTA 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S] 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

A  DISTRIBUTION  UNLIMITED:  PB  Public  Release 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

We  introduced  a  smooth  complex  for  nonlinear  elasticity  that  can  be  considered  as  the  tensorial  analogue 
of  the  standard  grad-curl-div  complex.  This  mathematical  structure  simultaneously  describes  the  kinematics 
and  the  kinetics  of  large  deformations.  The  relation  between  this  complex  and  the  de  Rham  complex  allows 
one  to  readily  derive  the  necessary  and  sufficient  conditions  for  the  compatibility  of  displacement  gradient 
and  the  existence  of  stress  functions  on  non-contractible  bodies.  The  main  application  of  the  nonlinear 
elasticity  complex  is  in  developing  mixed  finite  element  methods  for  large  deformations,  which  will  be 
pursued  in  a  future  project.  To  this  end,  the  smooth  complex  should  be  extended  to  also  include  less 
smooth  tensors.  We  introduced  this  extension  by  using  the  so-called  partly  Sobolev  spaces.  The  result  is  a 

Hilbert  complex  involving  second-order  tensors  on  flat  compact  manifolds  with  boundary.  We  then  used  the 
general  framework  of  Hilbert  complexes  to  write  Hodge-type  and  Helmholtz-type  orthogonal 
decompositions  for  second-order  tensors.  As  some  applications  of  these  decompositions  in  continuum 
mechanics,  one  can  study  the  strain  compatibility  equations  of  nonlinear  elasticity  in  the  presence  of 

Dirichlet  boundary  conditions. 

15.  SUBJECT  TERMS 

geometric  theory  for  nonlinear  elasticity,  discrete  exterior  calculus 

U.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF 

Standard  Form  298  {Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.1£ 

DISTRIBUTION  A:  Distribution  approved  for  pubiic  reiease 


https://livelink.ebs. afrl.afmiFlivelink/llisapi. dll 


8/25/2015 


FORM  SF  298 


Page  2  of  2 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

UU 

Arash  Yavari 

Unciassified 

Unclassified 

Unclassified 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

404-894-2436 

standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std .  Z39 . 1 8 


DISTRIBUTION  A:  Distribution  approved  for  pubiic  reiease 


https://livelink.ebs. afrl.afmiFlivelink/llisapi. dll 


8/25/2015 


1  Summary 


1 


Final  Report  for 

GEOMETRIC  STRUCTURE-PRESERVING  DISCRETIZATION 

SCHEMES  FOR  ELASTICITY 

AFOSR  GRANT  NUMBER  FA9550-12-1-0290 
Arash  Yavari 

Department  of  Civil  and  Environmental  Engineering 
Georgia  Institute  of  Technology,  Atlanta 


1  Summary 

We  introduced  a  smooth  complex  for  nonlinear  elasticity  that  can  be  considered  as  the  tensorial 
analogue  of  the  standard  grad-curl-div  complex.  This  mathematical  structure  simultaneously  describes 
the  kinematics  and  the  kinetics  of  large  deformations.  The  relation  between  this  complex  and  the  de 
Rham  complex  allows  one  to  readily  derive  the  necessary  and  sufficient  conditions  for  the  compatibility 
of  displacement  gradient  and  the  existence  of  stress  functions  on  non-contractible  bodies.  The  main 
application  of  the  nonlinear  elasticity  complex  is  in  developing  mixed  finite  element  methods  for  large 
deformations,  which  will  be  pursued  in  a  future  project.  To  this  end,  the  smooth  complex  should 
be  extended  to  also  include  less  smooth  tensors.  We  introduced  this  extension  by  using  the  so-called 
partly  Sobolev  spaces.  The  result  is  a  Hilbert  complex  involving  second-order  tensors  on  flat  compact 
manifolds  with  boundary.  We  then  used  the  general  framework  of  Hilbert  complexes  to  write  Hodge- 
type  and  Helmholtz-type  orthogonal  decompositions  for  second-order  tensors.  As  some  applications 
of  these  decompositions  in  continuum  mechanics,  one  can  study  the  strain  compatibility  equations  of 
nonlinear  elasticity  in  the  presence  of  Dirichlet  boundary  conditions. 

For  developing  finite  element  methods  for  nonlinear  elasticity  using  the  above  Hilbert  complex,  we 
first  discretize  this  complex  by  means  of  appropriate  finite  element  spaces  and  drive  stable  mixed  finite 
element  methods  for  the  associated  Laplacian.  We  use  the  general  theory  for  approximation  of  Hilbert 
complexes  and  the  finite  element  exterior  calculus  and  introduce  some  stable  mixed  finite  element 
methods  for  the  Laplacian  of  second-order  tensors  with  appropriate  mixed  boundary  conditions.  To 
this  end,  we  introduce  new  finite  elements  for  second-order  tensors  that  are  the  tensorial  analogues 
of  some  standard  finite  elements  for  vector  fields.  One  important  feature  of  the  finite  element  spaces 
generated  by  these  finite  elements  is  that  they  respect  the  global  topology  of  the  underlying  domains  in 
the  sense  that  they  reproduce  some  topological  properties  of  the  domains  regardless  of  the  refinement 
level  of  meshes. 


2  Nonlinear  Elasticity  Complex 

Let  H  c  be  an  open  subset  and  suppose  {X^}  is  the  Cartesian  coordinates  on  B.  We  equip  B  with 
metric  G,  which  is  the  Euclidean  metric  of  M^.  The  gradient  of  vector  fields  and  the  curl  and  the 
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divergence  of  (Q)-tensors  are  defined  as 

grad  :  X(i3)  ^  (grad  =  Y\j, 

curl :  T{®^TB)  T{®^TB),  (curlT)^-^  =  eiklT'^^ ,k, 

dW-.Ti^^TB)  ^X{B),  (divT)^  = 

where  ‘(j”  indicates  dldX"^ .  We  also  define  the  operator 

curF  :  T{^^TB)  T{^^TB),  (curf  =  (curlT)-^^ 

It  is  straightforward  to  show  that  curl"''  o  grad  =  0,  and  divocurl"''  =  0.  Thus,  we  obtain  the  following 
complex 

0  ^  x(i3)®-^r((2)2ri3f^Y((2)2ri3)  ^  x{B)  o,  (2-i) 

that,  due  to  its  resemblance  with  the  gcd  complex  (i.e.  the  standard  grad-curl-div  complex  of  vector 
fields),  is  called  the  gcd  complex.  Interestingly,  similar  to  the  gcd  complex,  useful  properties  of  the 
gcd  complex  also  follow  from  the  de  Rham  complex.  This  can  be  described  via  the  M^-valued  de  Rham 
complex  as  follows.  Let  d  ■  VL^{B)  ->•  be  the  standard  exterior  derivative  given  by 

(d(3)lo-lk  = 

i=0 

where  the  hat  over  an  index  implies  the  elimination  of  that  index.  Any  a  e  M^)  can  be  considered 

as  Q  =  (a^,  Q^),  with  od  e  ^}^{B),  i  =  1,2,3.  One  can  define  the  exterior  derivative  d  :  17^(13; M^)  ^ 

^  doP' ^  Since  dod  -  0,  we  also  conclude  that  dod  =  0,  which  leads  to  the 

M^-valued  de  Rham  complex  (I7(R; M^),  d).  Given  a  e  17^(13; M^),  let  [o;]*/!--/*,  denote  the  components 
of  cd  e  Q^{B).  By  using  the  global  orthonormal  coordinate  system  {X^},  one  can  define  the  following 
isomorphisms 


*0 

X(i3)  n^{B;R^), 

[2o(l^)]'  = 

*1 

*2 

[*2(T)]*jA'  =  ^u^jklT^^ 

*3 

X{B)  ^  n^{B;R^), 

N(1^)]'i23  = 

where  8j/  is  the  Kronecker  delta.  Let  T"*"  be  the  transpose  of  T,  i.e.  and  let  {E/}  be 

the  standard  basis  of  M^.  For  T  e  T(0^TB),  we  define  Tn  to  be  the  traction  of  T"*"  in  the  direction  of 
unit  vector  N  =  A''^E/  e  S^,  where  c  is  the  unit  2-sphere.  Thus,  Tn  =  N^T^'^'Ej.  One  can  write 

,4(^E2)  ,  A:  =1,2.  (2.2) 

It  is  easy  to  show  that 

o  grad  =  d  o  ^o,  12°  curl^  =  d  o  i^o  div  =  d  o  *2- 
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Therefore,  the  following  diagram  commutes  for  the  gcd  complex. 


0 - ^  X{B)  T{^‘^TB)  X{B) - ^  0 


^Q 


*3 


0  ^  n°{B-,  M^)  A  n^{B;  M^)  A  n^{B;  M^)  A  n^{B-,  M^)  ^  0 


(2.3) 


The  contraction  {T,Y)  of  T  e  r(®^Ti3)  and  Y  e  X{B)  is  a  vector  field  that  in  the  orthonormal 
coordinate  system  {X^}  reads  {T,Y)  =  Clearly,  if  Nq  is  the  unit  outward  normal  vector 

field  of  a  closed  surface  C  c  B,  then  {T,Nc)  is  the  traction  of  T  on  C.  Suppose  is  the  A:-th 

cohomology  group  of  the  gcd  complex.  Diagram  (2.3)  implies  that  ik  also  induces  the  isomorphism 
Hgcd('^)  ~  between  the  cohomology  groups.  Using  this  fact,  one  can  prove  the  following 

theorem. 

Theorem  2.1.  An  arbitrary  tensor  T  e  T{®^TB)  is  the  gradient  of  a  vector  field  if  and  only  if 

curfT  =  0,  and  £{T,ti)dS  =  0,  ^icB,  (2.4) 

where  I  is  an  arbitrary  closed  curve  in  B  and  ti  is  the  unit  tangent  vector  field  along  £. 


Similarly,  one  can  derive  the  necessary  and  sufficient  conditions  for  the  existence  of  a  potential  for  T 
induced  by  curl"''.  The  upshot  is  the  following  theorem. 

Theorem  2.2.  Given  T  e  T{0^TB),  there  exists  W  e  T{0^TB)  such  that  T  =  curl'''VU,  if  and  only 

if 

divT  =  0,  and  £{T,Nc)dA  =  0,  VC  c  13,  (2.5) 

where  C  is  an  arbitrary  closed  surface  in  B  and  Nc  is  its  unit  outward  normal  vector  field. 


We  can  also  write  an  analogue  of  the  gcd  complex  for  two-point  tensors.  Let  5  =  with  coordinate 
system  {x*},  which  is  the  Cartesian  coordinates  of  M^.  Suppose  99  : 13  ->  5  is  a  smooth  mapping  and  let 
Tx^p{B)  ■=  T^(^x)S.  Note  that  although  ip  is  not  necessarily  an  embedding,  the  dimension  of  Txp{B) 
is  always  equal  to  dim 5.  We  can  define  the  following  operators  for  two-point  tensors  that  belong  to 
T{Tp{B))  and  T{Tp{B)  ®TBy. 

Grad  :  T{Tp{B))  V{Tp{B)®TB),  (GradD)*^  =  U\i, 

CurF  :  T{Tp{B)  ®  TB)  ^  T{Tp{B)  ®  TB),  {Curl^ Ff  =  eiKLF"^,K, 

Div  :  T{Tp{B)  ®  TB)  ^  T{Tp{B)),  (DivF)*  = 

We  have  Curl"''  o  Grad  =  0,  and  Div  o  Curl"'"  =  0.  Thus,  the  GCD  complex,  which  we  also  call  the 
nonlinear  elasticity  complex,  can  be  written  as; 

0  ^  r(r(^(i3))^'^r(r(^(i3)  ®  r(r(^(i3)  (8)  ri3)  ^  r(r^(6))  ^  0. 
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By  using  the  following  isomorphisms 

11  :  0  TB)  ^  n\B;M.^), 

12  ■■  T{T^{B)  0  TB)  ^  n^{B;R^), 
h:TiT^{B))^n\B-,R^), 

one  concludes  that  the  following  diagram  commutes. 


[io{u)Y  =  u\ 

[l2iF)YjK  =  £jklT^^, 
[I3{U)]\23  =  U\ 


O^T{Tip{B))^-^'^T{T<p{B)  ®TB)'^-^T{Tip{B)  ®TB)^T{T^{B))  ^0 


^Curf 


Div  1 


0^n^{B;R^)—^n\B-,R^) 


->  n‘^{B;  M^)  — ^  n^{B;  E^)  ^  0 


The  above  isomorphisms  also  induce  an  isomorphism  Hq(-q(13)  rs  where  Hgf.p(13)  is  the 

A;-th  cohomology  group  of  the  GCD  complex.  Let  {E/}  and  {ej}  be  two  copies  of  the  standard  basis 
of  M^.  For  F  €  r{Tip{B)  0  TB),  and  n  =  n*ej  e  §^,  let  Fn  =  iTF^^Fj  e  X{B).  Then,  one  can  write 

Ik{F)  =  (^tk(^Fei^  ,tk(^Fe2^  ,  k  =  1,2. 

Let  {F,Y)  ■■=  F^^Y^Si.  The  above  relations  for  the  GCD  complex  allow  us  to  obtain  the  following 
results  that  can  be  proved  similarly  to  Theorems  2.1  and  2.2. 

Theorem  2.3.  Given  F  e  T{Tip{B)  ®TB),  there  exists  U  e  T{T(p{B))  such  that  F  =  GradC/,  if  and 
only  if 

CurfF  =  0,  and  £{F,te)dS  =  0,  \f£cB. 

Moreover,  there  exists  'I'  e  T{Ttp{B)  0  TB)  such  that  F  =  Curl"''’®',  if  and  only  if 

DivE  =  0,  and  J {F,Nc)dA  =  0,  MCaB. 


3  Hilbert  Complexes  and  Orthogonal  Decompositions 

Consider  the  following  linear  subspaces  of  V{Ti{){B))  and  T{T{p{B)  0  TB): 

T,{T^{B))  :=  {u  €  r{Tip{B))  :  U\g^s  =  «}  , 
rn/r99(i3)0ri3)  :=  {f  e  T{Tip{B)®TB)  :  Fe,  ±  OjB,  i  =  l,...,n], 
Tt,{T^{B)®TB)  :=  {f  e  T{Tip{B)®TB)  :  FeJ|c»,F,  i=l,...,n]. 


DISTRIBUTION  A:  Distribution  approved  for  pubiic  reiease 


3  Hilbert  Complexes  and  Orthogonal  Decompositions 


5 


The  operators  Grad,  Curl"'',  and  Div  can  be  restricted  to  the  above  subspaces  which  allows  one  to 
impose  boundary  conditions  on  the  GCD  complex.  The  upshot  is  the  following  commutative  diagrams. 


-  Grad;  —  —  Curl  _  _  Div,  - 


*0 


dn 


dn 


0  ^  nliB;  r3)  nUB-,  m^)  q2(b-,  m^)  nim  m^)  ^  o 


-  Divo  -  -  Curl  _  _  Grad ;  - 

0^T{Tip{B))  ^Tt^{Tip{B)®TB)  ^^T^^{T^{B)®TB)  ^fj{T^{B))^0 


-to 


0  ^  (S;R3)  nl(B0)  e- 


-I3 


nUB-,  r3)  nUB;  r3)  ^  0 


The  isomorphisms /o)"') -^3  are  L^-isometries.  The  Hilbert  spaces  L^r(r<^(i3)),  H^T j{T ^p{B)) ,  ^-{T {p{B)Si 

TB),  and  H^rtj{Tip{B)e)TB)  are  the  completions  of  {T{T(p{B)),  {{,  ))l2),  {Tj{T(p{B)),  {{,  ))Hi),  {Bnj{Tip{B)'S> 
TB),Uho),  and  (rt/T<^(i3)0rH),((,  ))ho),  respectively.  These  Hilbert  spaces  allow  one  to  write 
the  following  Hilbert  complex  for  two-point  tensors: 


0  ^  H^ri{Tip{B)f^"H^n^{Tip{B)®TB)^-^ 

H^ti{Tip{B)®TB)^" L‘^T{Tip{B))  0 


(3.1) 


The  dual  of  this  Hilbert  complex  reads: 

0  ^  L‘^r{Tip{B))?^"H^t2{Tip{B)®TB) 


(3.2) 


The  complex  (3.1)  is  isomorphic  to  {H^Vla_^{B) ,  dj^^) ,  and  hence,  it  is  Fredholm  with  Hg^-p  IB)  . 
Hqcdi('^)  ~  ^1'^)’  where  HQ(-p^(i3)  and  HQ(-p^(i3)  are  the  A:-th  cohomologies  of  the  smooth 

GCD  complex  (with  boundary  conditions  on  diB)  and  the  Hilbert  complex  (3.1),  respectively.  Let 
■■=  kerCurl]"  nicer Div2  be  the  kernel  of  the  Laplacian  L,^  associated  to  (3.1)  and  (3.2). 
Then,  T-L^^  consists  of  smooth  harmonic  two-point  tensors  and  “  ^gcDi('^)- 

can  show  that: 

Theorem  3.1.  Let  B  c  be  a  smooth,  compact  3-manifold  with  boundary  and  suppose  y?  :  H  ->•  R^ 
is  a  smooth  mapping.  The  Hilbert  complex  (3.1)  induces  the  following  L'^ -orthogonal  decompositions: 
The  Hodge  decomposition 

L'^r{Tip{B)  0  TB)  = 

Grad{H^T,{Tg,{B)))  ©  n^JB)  ©  Curf  {H^Tr,,{Tg,{B)  ®TB)), 

and,  equivalently,  the  Helmholtz  decompositions 

L‘^r{Tip{B)  0  TB)  =  Grad{H^ri{Tip{B)))  ©  kerDivs 

=  ker  Curl]"  ©  Curl^(H'‘^rn2(rv?(H)  0  TB)) , 
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where 


kerDiv2  = 

kerCurl7  =  Gra.d{H^Ti{Tip{B))) 

If  in  addition  a  two-point  tensor  is  of  class  (C°° ),  then  its  components  in  the  above  decompositions 
are  of  class  C"’’^  (C°° )  as  well. 

Corollary  3.2.  Let  i3  c  be  a  smooth,  compact  3-manifold  with  boundary  and  suppose  ip  :  B  ^ 
is  a  smooth  mapping.  The  necessary  and  sufficient  conditions  for  the  existence  of  a  Gradi -potential 
for  F  €  L‘^T{T<p{B)  0  TB)  are 

F€H^r^fTip{B)®TB),  CurfF  =  0,  {{F,K))l2  =  0, 

Similarly,  the  necessary  and  sufficient  conditions  for  the  existence  of  a  CurlJ  -potential  for  F  are 

F€H^t,{Tp{B)®TB),  DivF  =  0,  {{F,K))l2  =  0,  'IK  ^n^^  jB). 

Similar  results  are  also  valid  for  2-manifolds  in 


4  Mixed  formulations  of  the  tensor  Laplacian 


By  the  (weak)  3D  Laplacian  for  second-order  tensors  (or  simply  the  3D  tensor  Laplacian)  with  mixed 
boundary  conditions,  we  mean  the  following  boundary- value  problem:  Given  Q  e  L^{^^TB),  find 
T  €  LP‘{®^TB)  such  that 


■^{T)  =  Q, 

T  ±  Si,  (divT)ls,  =  0,  T||S2,  (curf  T)  ±  S2, 


(4.1) 


where  .if  :=  curl"''  o  curl"''  -  grad  o  div.  The  space  of  solutions  of  (4.1)  for  Q  =  0,  is  the  space  of 
harmonic  tensors  defined  as 


H{B,§i,§2)  ■=  ker.ifi  =  kercurlj^  nkerdivs2. 


One  can  write 

dimn{B,§i,§2)  =  dim5{^(i3,Si)  =  3bi{B,§i).  (4.2) 

Note  that  H{B,§i,§2)  measures  the  non-uniqueness  of  solutions  of  (4.1)  in  the  sense  that  if  T  is  a 
solution  of  (4.1),  then  so  is  T  +  T,  IT  e  ?^(i3, Si, S2).  The  integer  bi{B,§i)  depends  on  the  topological 
properties  of  both  B  and  Si.  For  example,  if  is  a  3D  ball  with  a  spherical  hole,  then  bi{B,0)  =  0, 
and  bi{B,dB)  =  1.  Therefore,  (4.2)  provides  a  connection  between  solutions  of  (4.1)  and  topological 
properties  of  B  and  Si. 

Since  'H{B,$i,$2)  is  finite-dimensional,  it  is  a  closed  subspace  of  L^{<0^TB),  and  therefore,  using 
the  orthogonal  projection  theorem,  one  can  write  L^(0^TB)  =  'H{B,§i,§>2)  ®  'H{B,§i,§>2)'^,  where 
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is  the  orthogonal  complement  of  i-e. 

n{B,§i,S2y  ■■=  {T  e  L\®^TB)  :  ((T,  S))l2^2  =  0,  VS  e  ^(i3,Si,S2)}  • 

One  can  show  that: 

Theorem  4.1.  The  boundary-value  problem  (4-1)  admits  a  solution  if  and  only  if  Q  €  'H{B,$i,$2)^ ■ 


For  approximating  solutions  of  the  tensor  Laplacians  with  mixed  boundary  conditions  by  means 
of  Galerkin  methods,  we  need  suitable  weak  formulations  for  these  boundary-value  problems.  The 
boundary- value  problem  (4.1)  can  be  considered  as  a  special  case  of  the  abstract  Hodge  Laplacian 
and  therefore,  one  can  define  the  following  mixed  formulation  for  (4.1):  Given  Q  e  find 

iU,T,P)  e  H\TB,§i)  x  H^{®‘^TB,$i)  x  H{B,§>u%2)  such  that 

iU,r))L2T  -  ((T,gradT))i2^2  =  0,  VT  e  H\TB,$i), 

((gradf7,0))i2^2  +  ((curl'^T,  curl'^0))i2^2+((P,  ©))i2^2 

=  ((Q,0)W,  y@€H^{^^TB,§i), 

((T,n))i2^2  =  o,  vn€?^(s,Si,S2). 

The  above  equations  are  the  Euler-Lagrange  equations  corresponding  to  a  saddle  point  of  the  functional 
:  H^{TB,§i)  X  H^{®‘^TB,$i)  x  H{B,$i,$2)  ^  M  given  by 

^(T,0,n)  =  ^((T,T))i2r-  ((gradT,0))^2^2  -  ^((curf  0,  curf  0))^2^2 
-  ((n,  0))i;,2(g,2  +  ((Q,  @))l2^2. 

The  first  equation  in  (4.3)  implies  that  U  =  divgjT.  The  second  equation  says  that  T  solves  the 
equation  ^{T)  =  Q  -  P,  where  the  harmonic  tensor  P  is  the  harmonic  part  of  Q  provided  by 
a  Hodge-type  decomposition.  Hence  Q  -  P  €  H.{B,§i,§>2)'^ ,  which  is  the  necessary  and  sufficient 
condition  for  the  existence  of  a  solution  due  to  Theorem  4.1.  Finally,  the  third  equation  eliminates  the 
degree  of  freedom  for  choosing  a  solution  T  by  requiring  that  T  e  ^(H,  §i,  S2)'‘'.  The  well-posedness 
of  the  mixed  formulation  (4.3)  follows  from  the  Fredholm  property  of  the  associated  Hilbert  complex. 
The  upshot  is  the  following  theorem. 

Theorem  4.2.  The  mixed  formulation  (4-3)  is  well-posed.  Thus,  for  any  Q  e  L^{®^TB) ,  the  problem 
(4-3)  admits  a  unique  solution  {U,  T,  P)  e  H^{TB,  Si)  x  H^{®^TB,  Si)  x  'H{B,  Si,  S2),  and  there  is  a 
constant  c  >  0  such  that 

+  ll^lliy=  +  ll-P|lL2,gi2  <  c||(5||i^2,g|2. 


In  the  mixed  formulation  (4.3),  the  boundary  conditions  on  Si  and  S2  are  the  essential  and  the 
natural  boundary  conditions,  respectively,  i.e.  the  boundary  conditions  on  Si  are  explicitly  imposed 
in  the  solution  spaces  while  the  boundary  conditions  on  S2  are  imposed  by  the  mixed  formulation. 
Alternatively,  one  can  write  another  mixed  formulation  of  (4.1)  as  follows:  Given  Q  e  L^{^^TB),  find 
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{S,  T,  P)  e  §2)  x  H^{^^TB,  §2)  x  n{B,  Si,  §2)  such  that 

-  ((T,curf  S))^2^2  =  0,  Vi:  e  H^{^^TB,§2), 

((curFS,  0))i2^2  +  ((divT,div0))i^22.+((P,0))i;,2^2 

=  ((Q,©)W,  V0eP'i(®2T^,S2), 

((T,n))^2^2  =  0,  VneP(i3,Si,S2). 

In  this  mixed  formulation,  in  contrary  to  (4.3),  the  boundary  conditions  on  S2  are  the  essential 
boundary  conditions  and  those  on  Si  are  the  natural  boundary  conditions. 


5  Discrete  Hilbert  Complexes 

It  is  possible  to  discretize  the  Hilbert  complex  mentioned  earlier  using  some  appropriate  finite  element 
spaces.  Let  Bh  be  a  3D  triangulation.  Then,  one  can  write  the  following  polynomial  complexes  for 
second-order  tensors: 

0^V^{TBh,Si)^V:_^{^^TBh,§i)^Vt2{^^TBh,Si)^Vr-3{TBh)^0,  (5.1a) 

0^v!:^{TBh,§i)^V^_A®^TBh,§i)^^VtA^^TBh,§i)^Vr-2{TBh)^0,  (5.1b) 

{TBh,§i)^Vr{^‘^TBh,§i)  ^^Vti{^^TBh,§i)^Vr-2{TBh)^0,  (5.1c) 

{TBh,§i)^Vr{^‘^TBh,§i)  ^^Vti^^TBh,§i)  ^Pr-i(m)^0.  (5.1d) 

Note  that  (5.1a)  is  a  subcomplex  of  (5.1b),  (5.1b)  is  a  subcomplex  of  (5.1c),  and  (5.1c)  is  a  subcomplex 
of  (5. Id).  These  discrete  complexes  are  Hilbert  complexes.  Regardless  of  the  refinement  level  of  the 
mesh  Bh,  the  cohomology  groups  of  the  above  discrete  Hilbert  complexes  are  the  same  as  those  of  the 
nonlinear  elasticity  complex. 


6  Stable  Mixed  Finite  Element  Methods  for  the  Tensor  Laplacian 


The  bounded  cochain  projections  between  the  nonlinear  elasticity  complex  and  the  discrete  complexes 
introduced  earlier  allow  one  to  use  the  general  theory  for  approximation  of  Hilbert  complexes  for 
approximating  the  nonlinear  elasticity  complexes.  Consider  the  following  mixed  formulation,  which  is 
associated  to  the  first  Laplacians  of  the  discrete  Hilbert  complexes  (5.1):  Given  Q  e  Lp‘{®^TB),  find 
{Uh,  Th,  Ph)  e  Vh  X  Vh  ^  '^h  such  that 


{{Uh,  r))L2T  -  {{Th,  grad  T))i2^2  =  0,  VT  e 

((gradl7fe,0))i2^2  +  ((cur^T/i,  curF0))i2^2  -h  ((P^,  0))i2^2  =  {{Q,&))l2^2,  V0  e  V^,  (6.1) 

{{Th,'n.))L2^2  =  0,  VlleP/i, 


where  the  pair  ,  V^)  can  be  either 


{V^^^{TBh,2>i),V^A®^TBh,2>i))  or  {TBh,^i),Vr {®^TBh,^i)) , 
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with  r  >  1.  Also  depending  on  the  choice  of  {Yh  '^h  •=  kercurlj^  nkergrad)^,  is  the  space  of 

discrete  harmonic  tensors  associated  to  one  of  the  discrete  Hilbert  complexes  (5.1).  Clearly,  we  have 
c  H^{TB,§>i)  and  c  However,  T-Lh  i-  TliB, §>iY2),  in  general.  Therefore,  the 

discrete  mixed  formulation  (6.1)  is  a  generalized  (or  non-conformal)  Galerkin  method  for  the  mixed 
formulation  (4.3). 

The  Fredholm  property  of  the  original  Hilbert  complex  and  the  existence  of  the  bounded  cochain  pro¬ 
jections  imply  that  for  a  shape  regular  family  of  triangulations  {Bh},  the  mixed  finite  element  methods 
based  on  (6.1)  are  stable  and  convergent.  The  rate  of  convergence  is  determined  by  the  smoothness  of 
the  data  Q,  the  smoothness  of  the  domain,  and  the  degrees  of  polynomials  that  generate  the  finite  ele¬ 
ment  spaces.  One  can  also  show  that  the  error  is  the  optimal  order  allowed  by  discrete  solution  spaces 
if  there  is  sufficient  elliptic  regularity.  For  example,  let  (Vf^Yh)  -  ■, 

and  suppose  that  the  solutions  T  and  P  are  of  Sobolev  class  .  Then,  one  can  write 

11^  “  +  h\\grad{U  -  U /j)  11^2^2  +  h\\T 

+  Y\\curl^{T-Th)\\L^^2  +  h\\P-  Ph\\L2^2  =  0{Y^^). 

Therefore,  regarding  the  degree  of  the  polynomial  approximation,  all  components  converge  with  the 
optimal  order. 

For  implementing  (6.1),  one  needs  to  calculate  the  space  of  discrete  harmonic  tensors  Ph-  To  this  end, 
consider  the  following  homogeneous  problem:  Find  {Uh,Th)  e  x  such  that 

{{Uh,  r))L2T  -  {{Th,  grad  Y))i2^2  =  0,  VT  e  , 

{{gradUh,&))L2^2  +  ((curl'^T/j,curl'^0))i2^2  =  0,  V0  e 

Then,  {Uh,Th)  is  a  solution  if  and  only  if  =  0,  and  T/j  e  71^. 
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